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Abstract 


A numerical simulation program is being coilstrubted having the folldwirtg 
features: (1) infinite circular cylindrical geometry with angle -dependence, (2) in- 
clusion of incident particles, photoelectrons, secondary electrons, backscattered 
electrons, any gun emissions, and any internal current pathways including surface 
conductive layers, (3) ’^quasistatlc time-dependent iteration”, in which sheath . 
potential changes during particle transit times are ignored, (4) use of afiprokimate, 
locally-dependenl space charge density expressions in solving Poisson's equation 
for sheath potentials, with use of nufneriCa" orbit -following to determine surface 
currents, (5) incident particle velocity distributions isotropic or beam-like, or 
some superposition of these, flatlonales for each of these features are discussed. 


1. I^TROUMCTIOS 

The asymmetry between sunlit and shaded areas of a synchronous spacecraft 
Is a key feature of the differential spacecraft charging problem at synchronous 
altitude. A realistic numerical model for the plasma shOath surrounding a 


/ 


synchronous sjjocecraft must therefor© bO at least two-dlmenstonnl. The only 
ekistln^ two-dimensional simulation which Is completely sclf-conslstonl is that of 
Soop, ^ who did a tlme-depondent treatment for a sphere, in which several thousand 
photoeiectrons were followed numerically. Such time -dependent treatments have 
until nmi# provided relatively low accuracy for a given computational expense, al- 
though there now exist Improved interpolation techniques for deducing space charge 
and fluX from a limited amount of orbit Information, which may change this situation 
in the future. 

Tivo other more simplified treatments are noteworthy. SChr'oder^ assumed 
that photoelectron emission was spherically symmetric, and thereby obtained a 
self-conslstent solution for a urtlpotentlal sphere, which showed the presence of 
potential minima due to photoClectron space charge in some circumstances. 

Lafon^ assumed spherical or cylindrical symmetry for space charge due to ambient 
particles, and negligible perturbation of this symmetry by photoeiectrons. He thus 
obtained radially symmetric self-consistent sheath potentials, but angle -dependent 
photoelectron density profiles, again for unlpotentlal spheres and cylinders. 

Here we describe a two-dimensional self-consistent simulation which avoids a 
completely time -dependent treatment, but instead is based oh a "quasistatlC time- 
dependent" iteration described in Section 2.2. Althbu^ results from three-dimen- 
sional simulations are likely to become available In the near future, ^ it Is generally 
true that *^h6 simplest realistic simulations are advantageous in elucidating basic 
physical effects, whereas more complicated ones are most useful for quantitatively 
predicting detailed interactions. 


2. FEATURES OF NUMERICAL SPACECRAFT - CHARGING MODEL 

2.1 Infinite Circular Cylindrical Geometry with Angle-Dependence 

This geometry implies the use of a polar coordinate grid for computations. 
Several reasons for such a choice, in preference to the more obvious spherical 
geometry, are: 

(1) Although a spherical geometry, with rotational (azimuthal) symmetry 
about the 3 paCecraft- 3 Un axis. Is two-dlmertslottal in position space, it Is three- 
dimensional In velocity space because particles with different azimuthal angular 
momenta must be treated separately. 

(2) Many spacecraft are finite circular cyllttdere. 

(3) In spherical geometry with azimuthal symmetry, focusing of particles onto 
the spacecraft -sun axis occurs Ih some models, leading to singularities In fluxes 
and densities along this axis. SUch effects must be regarded as spurious since 
real spacecraft are uhltkely to have the high degree of symmetry necessary to 
produce them. 
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(4) An infinite eyiinder, having a aurfaee senior with distinct properties, can 
bo rotated with respect to the sunward direction to stua/ the effects of such rota- 
tion. in a spherical geometry with azimuthal symmetry, the corresponding surface 
feature would be an annulus about the spacecraft -sun axis, and no such rotation 
would be possible without destroying azimuthal symmetry, 

(5) one major feature of spherical as opposed to cylindrical geometry, that 
is, the more rapid decrease of potential with increasing radius, can be modeled 
in an approximate way by simply adding the appropriate fictitious contribution to 
^ ^ in Poisson's equation. 

2.2 Physical Prccesses 

The model la to include velocity distributions of; incident particles, photo- 
electroiis, secondary electrons, backscattCred electrons, and any gun emissions. 
Internal current pathways including surface conductive layers are also being 
included. 

2.3 Quosistaiic Time-Dcpcadent Iteraticin 

In this procedure, sheath potential Changes during particle transit times are 
Ignored. This leads to the following iteration scheme; An angle -dependent sur- 
face potential IS chosen. Poisson's equation Is then solved to provide a radius- 
and angle -dependent static sheath potential (see Section 2. 4 below). Particle 
orbits are then followed numerically ih this potential, yielding surface charging 
rate as a function of angle (orbit -following is, however, not used to provide space 
charge densities fbr Poisson's equation; see Section 2.4 below). These rates are 
then averaged over any conducting sector, and any currents transferred internally 
(including those through any surface conductive layers) are subtracted. The re- 
sulting net charging rates are then used to carry forward onetime step, yielding 
new surface potentials. This process is then repeated until a steady-state condi- 
tion results, or, in a situatloft Ift which external conditions vary with time, is 
repented to follow such time-dependence. 

■fhe use of this procedure, as opposed to a Completely time-dependent slmU- 
Iptton, should produce important computational economies. Clearly one will lose 
information about very rapid transient phenomena with thlts approach. However, 
steady-state or slowly time-varying situations are of major Importance, these 
include changes in the incident particle distributions, which are likely to have time 
scales of seconds or minutes. 
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2.4 o( ApproitM&te Spli«<*>ClMrge 04^rt^ity Expffif^Aicins 

At synchronous altitude, the Debye length Xq for ambient particles is usually 
> 10 m, so for satellites of Ordinary siae, effects of ambient space charge on 
sheath potentials will be relatively small. Any reasonably realistic approximation 
of this space charge can therefore be expected to produce only negligible errors in 
solving Poisson's equation for sheath potentials. Furthermore, large savings in 
computer time can be expected to resu:t if one can avoid exact density calculations 
involving numerical orbit -following. Ih the present Work, it is intended that a 
relatively Small amount of orbit-followihg be done to calculate surface currents 
(section 2.3). 

A more significant space-charge effect near the spacecraft may be caused by 
emitted photoelectrons or Secondary electrons, ^ because of their relatively low 
velocities compared to ambient values. However, effects of these are likely to 
also be small enough that any reasonably realistic approximations for their densi- 
ties wlli yield good accuracy. ^ Such approximations must ultimately be validated 
by comparison with a few carefully chosen exact calculations. It is advantageoxis 
if such approximations depend On local potential only (rather than potentials at 
many locations), together with a relatively small number of other parameters, 
sUch aS spacecraft potentials and potential barrier heights and locations. Here We 
propose three types Of space-charge density approximation, as follows. 

2. 4. 1 apphqximAtions foh poTentia l Wells without 
ohstacles 

Exact density expressions have been developed for CollisiohleSS, Maxwellian 
particles in the presence of obstacle-free potential wells of arbitrary shape by 
Laframboise and Parker. ^ The appropriate expression for our purposes is the 
result given by their Eq. (2) for three-dimensional wells. This is true even for 
art "infinite", that is, very long cylindrical Spacecraft geometry, because of par- 
ticle entry at the ends of such a geometry. For definiteness, We consider a 
negative well given by ^ (x, y, z) ^ 0, with 0 - 0 as + y^ + - oe , where ^ is 

electric potential. If only ambient particles are considered, PoissOn's equation is: 

2^=.^(N^.Nj) (1) 

O 

where e is magnitude of unit electron charge, is pormlttlvlty of space, and Ng, 
N{ are electron and ion number densities, respectively. Since positive ionS are 
the attracted species in this well, we use Eq. (2) of l<afranioolse and Parker^ for 
ion density, and the usual Holtzmcuin factor for electron density. If 
^De “ is electron or ioh density far from the spacecraft. 
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L Is & 6h&r&6tei^i3tic st>&ce6raft length, ^ = L^, x - c^/hTg< 0, k Is SoU^nlann'S 
constant and T> 14 xet»perattihe,,.SQ<...(Jl).he(iniueai . 

[l-xW‘'^ + e (2) 

Where g(s> = exp (s^)erfc(s) = exp(s^) j exp(-t^)dt. 

The important feature of Eq. (2) for oUr purposes is that its right-hand side 
is a function of x only. For small x, Eq. (2) reduces to; 

= (1 +tg/T.)(L/Xj>g)^X (3) 

Where terms of order x and higher have been ignored. The linear form of (3) 
permits the use of direct PoiSson-Solvers for finding x . Another simplified form 
ban be obtained by rederiviag Eq. (2) with monoenergetic instead of Maxwellian 
ions assumed. The appropriate monoenergetic velocity distribution {Chen;^ 
Laframboise , p. 1 4) is: 


__ 6(E - E^) 

d^v (2m.Ej)V2 


(4) 


where Ej^ =4kT./ff and m^ is ion mass; this distribution duplicates the ambient 
number density and flux values of a Maxwellian at temperature T.. Rederivation 
of (2) using this distribution yields the computationally simpler form: 



If any regions estlst where x > i). the roles of ions and electrons are inter- 
changed, and Etis. (2)-(S) must be modified accordingly. 

The essential approximation contained in Eqs. (2) - (5) is the neglect of orbit 
depletion due to intersection with the spacecraft. The densities of ambient ions 
and electrons will therefore both be overestimated near the spacecraft in these 
results. As long aS the spacecraft is at least moderately smaller than the 
effects of this overestimate will be small. The attracted-speCieS density Will be 
overestimated by the greater amo\mt for reasons involving the curvatures of 
attracted ind repelled particle orbits. The sheath profiles predicted by i2)l or 
(^) will therefore be steeper than real profiles, if electron emission effects are 
ignored. 
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2. 4. 2 APPROXIMATIONS EASED ON SVMMDtRIC. POTENTIALS . 

Laframbois^f, aftd Lafrarftbdlaa arid C3odA^*d, ® Eqa, (7) and (fi), have presented 
expressions for number densities of ambient attracted and repelled Maxvt^ellian 
particles, respectively, which ar*? exact for radially symmetr ic monotonlc poten- 
tials near a perfectly absorbing spherical collector. These expressions Contain 
terms identical to the ion and electron derlsity expressions in (2), together with 
subtractive terms representing the effects of particle interception by the Collector, 
Whipple® has used a thick-sheath approximation to develop density expressions for 
both ambient and emitted particles in the presence of a potential barrier, again for 
spherical symmetry. Lafon^ has developed approximate density expressions for 
escaping photoelectrons, based on assumed spherical or cylindrical symmetry in 
the sheath potential, but not in the photoemission fluxes. Since all of these expres- 
sions depend only on local potential and a small number of other parameters, it iS 
tempting to explore the possibility of using them even in the presence of sheath 
potentials which are known to be angle -dependent, and near spacecraft having nOn- 
spherlcal shapes. All of these expressions depend essentially on the solid angles 
subtended at any given radius, by orbits which have intersected the spacecraft, 
for all significantly populated particle energies, including ' • the effects of Orbit 
curvature due to electric fields. It is likely that in mahy caSes, such Solid angles 
will not be greatly modified by angular asymmetries in sheath potentials (from 
symmetry, such modification must be of second order in angular variations). In 
using such approximations With irregular spacecraft Shapes, It would be necessary 
to define some way of Choosing "radius” for substitution into them. One way to do 
this would involve matching the solid angle subtended by the spacecraft at the 
location in question, with that subtended by a sphere as a function of radius. 

Similar procedures would be necessary for dealing with parameters describing 
potential barriers irt these expressions, Lafon^®' and Parker^^ have given 
useful general discussions of the formulation of density expressions for symmetric 
potentials. 

2. 4. 3 APPROXIMATIONS BASED ON ECjUlVALENT POTENTIAL 

Wells 

We consider the idealized situation shown in Figure 1, in whi' h a spacecraft 
is assumed to have shaded-side surface potentials which are Very negative, and 
suhllt-side surface potentials Which are close to space potential. The solid curves 
outside the spacecraft represent ecjulpotentials. The dotted curve FGH represents 
a surface which passes through the saddle point G oii the sunlit side, and is every- 
where perpendicular to the equipoteritials, so this sUrface represents the maximum 
extent of a sunlit-side potential barrier for electrons. Fahleson^^ has pointed out 
that such a barrier may exist even wheh space charge is negilgible, because of the 
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Figure 1. G6ii6rsil Appearauiee pf a Possible 
Sheath Potehtlal Profile atPuhd a Spacecraft. 
Dotted curves inSlde the spacecraft surface are 
fictitiPUS extetisiohS of eqtilpotential surfaces 
outside, as described following Eq. (6) 




sunlit -Shaded asymmetry til surface potentislS. We consider aS en exSniple the 
process of approximately Calcillatlhg phPtPeleCtron space chsrge density inside 
this barrier; CSlculStion cf secondary electron charge density is similar In mpst 
respects. We consider all those phptoelectrohs emitted With a total energy Eg 
equal (within some differential amount dE) to the potential of the equipotential 
surface flBC. Siich particles can never go outside ASC, but must reimpihge on 
the spacecraft surface ADd. If lp(E) is the photOemisSion coefficient, that is, 
the energy-differential particle current density of photOemission from the space- 
craft surface (this will depend on surface material and solar illumination angle), 
then the totai-photoemissiOn particle current between energies Eg and EJg -i- dE is: 

s] (6) 

Where Fg + e^g > 0, s represents surface position, s ^(£g) IS surface potential, 
Fg + e^^ = Fg + e^^ = (), and the integration is over the surface ADC. Since t 
for most materials is largest for emission kinetic energies ^ higV^ = Eg+C(^y« 1 volt, 
most Of the photoemission between energies Fg and Eg + dF will tend to Come 
from regions such as. Say, A' and C In Figure 1, Where is about i volt more 
positive than at A and C. bn the other hand, particle motions Will tend to spread 
the reimpingeifient current more uniformly over ADC. 
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We now model this process epproxtm&tsly by mentally removing the spacecraft 
surface between A and C, and replacing It by an arbitrary extension AEC of the 
ec^nipotentlal surface ABC. We also do the same for other et|ulpOtentlals which lie 
Inside this one, as also shown in Figure 1. We have nOw" constructed" an obstacle- 
free potential well, and We can Us6 the Laframbolse-Farker® theory to derive 
model density and flux profiles for such a well. We cah then Integrate the latter 
over ADC and match the result with Bq. f6). We rewrite the mortoenergetiC 
distribution (4) for electrons as follows : 


m^N* 6(E-Eg) 


(7) 


N* is ho^ A reference ntiinbei' density to be evaluated. We obtain: 

N s Pf [ ^(r)] = y*fd^ V = N* [ 1 + e^ir)/Eg] H^eilr) + Eg] __ (8) 

J a J[ ^(r)] = f fV^ d^v = N*(Eg/8mg)^/^[ 1 + e(|t(r)/Eg] H[ e<^(r) + Eg] (9) 

Where j is a number flux erossing ah arbitrarllj^ oriented surface element from 

either direction, v is Velocity component jjerpendicuiar to such a surface element, 

/s ^ 

d(x) dx is the Heaviside step function. The total number flux 

-go 

Crossing ADC from either direction iS ftoW given for our model well by: 

. ( 10 ) 


our procedure for approximating the Space charge density now involves per- 
forming the integrations over the surface ADC in both (6) and- (10), then evaluating 
N* by equating these two results. 'This Is done for each of the discrete energies 
Egj which are chosen to represent the-photoemisslon. The quantity dE in (6) must 
then be chosen equal to the separation between these energies. The resulting set 
of values is then used together with (6) tb construct the space-charge density 
expression; 


N(^(r)) = t N*t l + e*(r)/Egjl^/2 H(e^(r) + Eg^] 
3 


(ID 


This expression also has the advantage of dependence bnly on local potential, 
as db those derived ih Sectlbhs 2.4. 1 and 2.4. 2. in Using it, one would precal- 
culate the coefficients N* as described above, then use (11) as a contribution to 
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the sp£ie6 charge density in Poisson's ei)uatlon. Ah important aptiroxiihatioA con- 
tained in (11) Involves negleot of the fact that photoemission fluxes given by (lO) 
are in general distributed differently o\ter ADC than those given by (6). For 
energies Ej > -etJ^ > 0, where is the saddle-point potential in Figure 1, some 
photoelectrons would escape, and the corresponding terms in (11) Would be over- 
estimates. 

2.5 USe-ol isotropic or Beam-Like Incident Velocity Distributions 

Important computational economies clearly result from assuming that incident 
velocity distributions are either isotropic or beam-like (monokinetic); the approxi- 
mate density expressions described in. Section 2. 4 are examples of results for 
isotropic distributions. Any incident distribution may be modeled as closely as 
desired by a superposition of isotropic and beam-like distributions. 


3. CONCLUSIONS 

We have described the major features of a "quaslStatlc time-dependent" 
numerical simulation of differential spacecraft Charging at synchronous altitude, 
incorporating an infinite cylindrical geometry v?lth angle -dependence. Ifhe com- 
puter program involved is presently under construction. 
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